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SUMMARY 

An electrolytic tank analog used for two-dimensional analyses of electrostatic - 
thrustor ion optics is described in detail. This analog has unique features that increase 
the ease of operation and versatility beyond previously reported designs. A generalized 
procedure, which permits a wide variety of space-charge-flow problems to be solved 
quickly and accurately, is used for space-charge simulation. 

Solutions are given for operational problems encountered in the practical use of the 
analog, and some analyses of specific thrustor accelerator configurations are presented 
to illustrate procedures and important applications of the analog. The use of the analog 
as a precursor to a more accurate digital computer analysis is recommended to speed 
up the initial investigation of a thrustor configuration. 


INTRODUCTION 


SaTf+ciZ 



Many of the problems associated with electrostatic thrustors have been solved. The 
focus of attention is turning toward lifetime extension (ref. 1). To maximize the opera- 
ting lifetime of a thrustor, the electrode system must be designed so that the ion beam 
does not intercept the accelerator electrode (see fig. 1). 

The width and shape of the ion beam in a thrustor is 
a result of the combined effects of the space charge and 
the electrode shapes. Analyzing the space-charge flow 
for a wide variety of thrustor configurations can be rel- 
atively simple and fast if an analog is available that is 
sufficiently versatile, accurate, and convenient to use. 
Digital computer programs exist for analyzing space- 
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Figure 1 . - Cross-sectional diagram of electro- 
static thrustor. 














charge flow in considerable detail, but an electrolytic tank analog is more appropriate 
for initial investigations of thrustor configurations. Ions produced by charge exchange 
can also be treated readily with the electrolytic tank so that it is possible to estimate 
electrode erosion and lifetime (ref. 2). 

The electrolytic tank analog discussed in this report is based on a design procured 
by Lewis from Hughes Research Laboratory. As will be shown, however, the equipment 
and techniques reported herein have been developed and extended to permit the solution 
of a wider variety of space-charge-flow problems in a shorter time. Along with a dis- 
cussion of these solutions and the methods used to obtain them, this report describes the 
basic considerations involved in the design and operation of the analog. Sample calcula- 
tions used in the operation of the analog are provided. 


APPARATUS AND PROCEDURE 


The plans and procedures for the original electrolytic tank analog at Lewis Research 
Center were procured from Hughes Research Laboratory. The present Lewis Research 
Center analog consists of the electrolytic tank, an analog computer, and a large X-Y re- 
corder for obtaining permanent records of the results (fig. 2(a)). This analog has, since 
procurement, been considerably modified in both equipment and procedures. A basic 
electrical diagram of the system is shown in figure 2(b). 

To simulate the thrustor for the analog investigation, scale-model electrodes and an 
ion-emitting surface are placed in a plastic tray filled with electrolyte. For symmetri- 
cal designs, only half the thrustor configuration need be considered. Voltages applied to 
the model electrodes establish electric fields throughout the electrolyte, and the space 
charge due to ion flow is simulated by means of "current” injection pins located in the 
tray (refs. 3 and 4). 

The space -charge-flow problem involves, of course, the solution of Poisson’s equa- 
tion (v V = p/e ) subject to the initial and boundary conditions of the particular case. 
(Symbols are defined in appendix A. ) The equations of motion for the two-dimensional 
problems are 


and 


d 2 x r 

"dt 2 ” 


d 2 y T 


dt 




= -t)E 


x,r 


= -rjE 


y>r 


( 1 ) 


3 



where 77 = q/m is the charge-to-mass ratio of the 

ions, E and E are the x- and y-components of 
X, a y, r 

the electric field intensity, and the subscript r refers 
to the real thrustor. The equations of motion for a 
similar particle in the analog model are obtained by 
dimensional and time scaling so that, for the model 
thrustor, 
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Figure 3. - Field-sampling probe details. 
(All dimension in inches.) 


where x m and y m are the scaled model dimensions, 
r is the scaled time, and E x m and E^ m are the 
scaled electric field components set up in the electro- 
lyte by the model. Equations (1) and (2) are of the 
same form and may be treated in the same manner. 

To begin the solution, the space-charge-free electric 
fields are used as the first approximation to the solu- 
tion of the space-charge-flow problem, and ion trajec- 
tories are determined under these conditions. These 
trajectories provide a means of estimating a charge 
distribution from which a better approximation for the 
trajectories may be determined. By continuing this 
iterative process, convergence to the solution of 
Poisson’s equation is obtained. In actual operation, 
the x- and y-components of the field intensity, meas- 
ured at the surface of the electrolyte by a set of sam- 
pling probes, are fed to an analog computer for the 
solution of the equations of motion. The sampling 
probe carriage, driven by the computer output, traces 
out ion trajectories. The X-Y recorder shown in fig- 
ure 2 is synchronized with the probe carriage and pro- 
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Figure 4. - Divergent-flow model electrodes in tank 
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vides a permanent record of the trajectories. Figure 3(a) 
is a photograph of the underside of the field- sampling probe. 
The probe details will be discussed in the section Field- 
Sampling Probe. 


Electrolyte Tray and Current-Injection Pins 


F Trnensions mTnch«.T tl0n P ' n ' (A " The P lastic electrolyte tray is approximately 16 inches 

wide, 30 inches long, and 2 inches deep. The floor of the 
tray is accurately leveled by adjustments at its supports. Small sockets sealed in the 
floor of the tray provide receptacles for removable current- injection pins (see figs. 4 
and 5). This unique feature permits the pins to be plugged in or removed as needed. 
There are several advantages of this removable-pin design compared with a permanent- 
pin design. The model electrodes do not require holes in them to fit over the pins. 
Cleaning the floor of the tray is much easier, and most important, electrode configura- 
tions may be changed quickly. Although trays exist in which the pins are secured by 
screws from the underside of the tray floor, the plug-in pins are more convenient to use. 
A patch-panel system (fig. 6) is used to connect each current- injection pin to a potenti- 
ometer and a series resistor. By adjusting the potentiometer and monitoring the voltage 
dropped across the resistor, the injection current of each pin may be set. The calcula- 
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Figure 6. - Current-injection patch panel. 


tion of the value of this current will be discussed in the section Obtaining Working Equa- 
tions for Analog. Samples of the calculation are provided in appendix B. 


Field-Sampling Probe 

A photograph of the field- sampling probe is shown in figure 3(a) (p. 4). The probe 
head is attached to a micrometer screw so that the probe may be lowered into the elec- 
trolyte to a precise depth. The probe head consists of a nylon mounting block with four 
small gold-plated brass pins protruding from the underside. The brass sampling pins 
are held in locating holes in the nylon block by set screws. Electrical connections to the 
pins are made through the set screws. The pins form a square array (fig. 3(b)) with the 
diagonally opposite pins spaced 1/10 inch apart. One of the brass pins is shown in detail 
in figure 3(c). References 5 to 7 discuss probe spacing and the effects of meniscus on 
the accuracy of electrolytic tank analogs. In principle, the closer the pin spacing is, the 
more accurate the field sampling will be; however, the distortion due to the meniscus 
about the pins may overshadow the increased accuracy due to closer spacing. 
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Associated Circuitry and Components 


A power amplifier, driven by a 100-cycle-per-second sine- wave generator, delivers 
power to a transformer (fig. 2(b), p. 2). The model electrode potentials and the bias 
voltages of the current- injection pins are obtained from the multitap secondary of this 
transformer. Alternating-current voltages rather than direct-current voltages are used 
to reduce the battery action or polarization that occurs between electrodes. 

The cathode followers (fig. 2) have a high input impedance to minimize current flow 
through the field- sampling probe pins; thus, field distortion by the pins is prevented. 

The output of the x- cathode followers is fed through a Wagner grounding circuit (fig. 2(b) 
and ref. 8). The output of the Wagner ground amplifier is then fed back to the model 
electrodes. As the probes move through the electrolyte, the output of the Wagner ground 
amplifier varies in such a way so as to maintain an electrical ground nearly midway be- 
tween the two x-probes, thereby reducing capacitive effects at the probes. This results 
in a shifting in the values of the electrode reference potentials; however, the net poten- 
tial difference between them remains unchanged. The values of the amplifier input ca- 
pacitor and the amplifier gain of the Wagner ground circuit are adjusted so that the am- 
plifier will remain stable at as high a gain as possible. For the analog reported herein, 
the maximum stable gain of the amplifier is about 1000. 

The outputs of the x-cathode followers are also fed to the analog computer where 
their difference is determined and multiplied by an alternating- current reference signal, 
and the result is filtered to yield a direct-current signal. The magnitude of this signal 
is proportional to the field intensity E . The same operations of finding the difference, 

A 

multiplying by a reference signal, and filtering are performed on the outputs of the 
y-cathode followers, and the two components of the equations of motion are solved simul- 
taneously by the computer. 

The circuit used for obtaining the difference, multiplying by the reference signal, 
and filtering (see fig. 7) eliminates the accessory equipment used in other analogs for 
this purpose (e. g. , ref. 7). Computer components are used exclusively. The operation 
is as follows: for the signals from each pair of probes, the difference is obtained by in- 
verting the sign (shifting the phase by 180°) of one and adding that result to the other sig- 
nal. The resulting signal is a sinusoidal function A ^ sin cot where A ^ is the ampli- 
tude, co is the frequency, and t is the time. By using an electronic multiplier, a ref- 
erence signal is then multiplied by this function. The reference signal may be of the 
form Ag(cot + 0 ) where 6 is the difference in phase between the working and the refer- 
ence signals. The resultant product is then 


A.An A. An 

[Aj sin cot] [A 2 sin(cot + e] = — - cos 6 cos (2cot + 6 ) 


( 3 ) 
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Figure 7. - Block diagram for analog computer program. 


The first term of the product is a constant for any given phase difference, while the sec- 
ond is a sinusoidal signal of twice the original frequency. The low-pass filter circuit 
(fig. 7) allows only the constant portion of the product to be passed. This direct-current 
signal is proportional to the electric field strength component. Simultaneous double in- 
tegration of these field components in the computer circuit yields the trajectory of a 
particle. The outputs of the double integration are used to drive simultaneously the 
x- and y-servosystems in the electrolytic tank, which propel the probe. 

If the equations of motion are electronically integrated simultaneously by the com- 
puter only once and the computer outputs are plotted one against the other, the results 
are electric field lines. By reversing the plotter inputs, lines normal to the electric 
fields (equipotential lines) are plotted (ref. 7). If the gradient throughout an area is ex- 
tremely small, the analog may not be able to plot the potential field accurately. In this 
case, the field must be plotted manually point by point. 
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Calibration 

Before use, it is necessary to calibrate the analog. As mentioned previously, the 
four field- sampling probes are arranged to form the corners of a square (fig. 3(b), p. 4). 
Perfect alinement of these probes is practically impossible. Small deviations can be 
compensated for, however, by an electrical correction made to the signal measured by 
one pair of probes. This is done by adding part of the signal from one pair to the other 
so that both signals are the same under identical field conditions. Of course, the great- 
est degree of accuracy is obtained when the correction signal is very small. It is there- 
fore necessary to aline the probes with care. In running an actual problem, the accuracy 
of the analog is also improved if the correction signal is taken from the pair of probes 
that are normal to the general direction of the trajectories. 


Electrolyte and Model Electrodes 

Two important factors that affect the amount of injection current are the conductivity 
ct of the electrolyte and the value of the fixed resistors (fig. 2, p. 2) in series with the 
injection pins. These two factors are independent of the model configuration and must be 
adjusted so that the desired range of currents can be injected with the power supply that 

o 

is used. In the present analog, ordinary tap water having a conductivity of from 2. 0x10" 
-2 

to 4. 0x10 mho per meter is used for the electrolyte. The fixed resistors have a re- 
sistance of 10 000 ohms. The electrode material used for the models is stainless steel. 
Details of tests that were made to evaluate various materials are presented in appen- 
dix C. 


Initial Conditions 

Because of physical limitations imposed by the size of the sampling probe, trajec- 
tories to be traced are started about 3/16 inch from the simulated emitter. Because of 
this limitation, the simulated charged particles must be given an initial velocity that 
corresponds to the velocity a real particle would have in the field at that location. This 
condition may be set by applying the initial velocity in terms of voltage, at the initial 
condition junction of the velocity integrators (see fig. 7). 
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Spacing of Current-Injection Pins 


The most accurate simulation of a space-charge field is accomplished by using the 
smallest injection-pin mesh spacing possible. A practical limitation, however, is the 
number of calculations and settings required to determine the pin currents from iteration 
to iteration. Because of the complex network of resistances between pins in the tank, 
some pin currents change as others are set. When a large number of pins are used, they 
may require resetting as many as three times for each iteration. The mesh size found 
suitable for a variety of problems was 0. 8 inch. The number of pins used in the prob- 
lems described in the section ILLUSTRATIVE EXAMPLES ranged from 25 to over 200. 


THEORY 

Space-Charge Simulation 

The relation between the electric field and the charge density in a thrustor is given 
by Poisson's equation in the form 



where E r is the electric field and is the charge density in the real thrustor. 


(4) 

Define 


where E m is the electric field in the electrolytic model, 
that 
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m m 
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where F is the scale factor from equation (6). The relation between the current density 
and electric field in the electrolyte model is given by 


i = aE m 
J m m 


( 8 ) 


where a is the conductivity of the electrolyte. By substitution, equation (7) becomes 
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By applying Gauss’ law and including the scaling factors of equation (6), equation (9) 
becomes 
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where Q f is the charge in a volume element in the real system, which must be simula- 

lated in the model system. The region 
through which the beams of each system 
pass may be divided into a number of these 
volume elements (fig. 8); the elements are 
parallelepiped in shape, and those in the 
model system, called simulation volumes, 
are scaled replicas of those in the real 
system. If the integrated current through 
the elementary volume is denoted by I, 
equation (10) can be written 
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Figure 8. - Cross section of thrustor showing cross sections of simu- 
lation volumes. 


This current must be injected into the 
simulation volume in the electrolytic-tank 
analog from an external source. Since 
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field measurements are taken at the sur- 
face of the electrolyte, the proper poten- 
tial distribution must also occur at this 
surface. It has been empirically deter- 
mined (ref. 3) that these conditions may be 
satisfied by injecting the current from a 
partially insulated pin protruding from the 
floor of the electrolyte tray. As shown in 
figure 9 the pins actually simulate the 
charge in the cylinders. If the cylinders 
are equal in cross-sectional area to the 
parallelepipeds, however, a good simulation is obtained at the surface of the electrolyte. 
The determination of the appropriate value of the injection current I of equation (11) 
requires a knowledge of Q r . The method used herein to calculate the values of Q r dif- 
fers considerably from previously reported methods and is a key factor in increasing the 
versatility of the Lewis Research Center analog. 

To begin, consider the real thrustor and a scaled thrustor the size of the electrolyte 
model. The relation between the currents emitted from similar elemental areas of the 
emitters in the two thrustors is 


/-Mesh area 
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Figure 9. - Simulation of charge by cylinders representing 
parallelepipeds. 


( 12 ) 

where is a constant and is the current in the scaled thrustor. The total current 
flow from the thrustors is divided into " stream tubes. M Each stream tube of current 
flow originates from one of the elemental areas on the emitting surface. To consider 
space-charge- limited current flow, the currents flowing in a stream tube of either 
thrustor may be approximated by employing a technique described in reference 9 in con- 
nection with a digital computer program for solution of space-charge-flow problems. 
This technique is referred to as the "plane diode approximation. " To determine these 
currents, an equipotentional line close to the emitting surface is selected as shown in 
figure 10, and it is assumed that it forms a parallel- plate plane diode of separation 
with respect to each elemental area of the emitting surface. By using the mean value 
of J| is calculated for each stream tube from Child's law for space-charge-limited 
flow 


J f - c 2 v, 3/2 a ‘ 


t 2 


(13) 
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where Cg = (4/9)e Q -^2 rj, Vj is the potential difference shown in figure 10, and 
elemental area. Similarly 


7 _ c xfifa 

J r c 2 v r 


From equations (13) and (14) it follows that 


V 3/2 
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and the proportionality constant of equation (12) is therefore 
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For each elementary simulation volume, summing over all the stream tubes in the 
volume gives 
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(16) 



(a) Trajectories passing through mesh area (top view of volume element). 



Q r= Z J r,i l i 
i 

where J r . is the current in the i 
stream tube and t^ is the interval of 
time (measured in the real system) that 

IL 

the ions of the 1 tube remain in the 
simulation volume. In actual practice, 
of course, the smaller the elements of 
emitting area and the smaller the sim- 
ulation volumes, the better the space- 
charge simulation will be. By using 
equation (15), equation (16) becomes 


Q 

^r 


y 3/2 

i^EVi 

V, /Z i 


(17) 


(b) Trajectory with superimposed 
timing marks. Number of time 
spaces in simulation volume, 4. 


Now let 0 = T./t. where r. is the 
11 1 ,th 


Figure 11. - Typical mesh area. 


scaled time that the ions of the i 
stream tube spend in a simulation vol- 
ume of the model- size thrustor. The 

value of 0 is determined from the scaling constants between the real and simulated 
thrustor and will be discussed in the next section. Then, by substitution, equation (17) 
becomes 


,3/2 


Qr- = 


P V f 3 / 2 0 i 


Z J f, i T i 


(18) 


In an actual case, the stream tubes are constructed so that the particle trajectories form 
the centerlines of the tubes (fig. 11(a)). Timing marks are superimposed on the trajec- 
tories at equal time intervals T. Thus, for any simulation volume, the number of time 
spaces in the simulation volume N y i may be counted for each stream tube (fig. 11(b)). 
Multiplying N . by the time per space (T) yields r. for each stream tube. When only 
a portion of a stream-tube section passes through a charge simulation volume, the con- 
tribution to the charge in that volume is obtained by prorating on an area basis. For 
example, in figure 11(a) stream tube 2 contributes entirely to the charge in mesh 2C 
(double cross-hatched area). The prorating factor Aj is 1. 0 for this case. Stream 
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tubes 1 and 3 contribute only partially to mesh 2C (single cross-hatched areas). Since 


T i = N v,i T 


using equation (18) results in equation (11) becoming 


1 = 


a F 2 V 3 / 2 T 




- — y j f iN ,a. 
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c Pi 


(19) 


Now a may be written as 


a = — — F 


( 20 ) 


By using this equation and the Child’s law relation for the current flowing in a stream tube 


J — ^ c i/2n V 3 / 2 * 

J f,i-- e oV 2 ^ V f ^2“ 
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equation (19) becomes 
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( 22 ) 


The emitting area a^ may be taken outside the summation by requiring that all stream 
tubes be the same width. This is the basic equation resulting from the new technique 
used to evaluate space charge. The currents in the stream tubes are recalculated for 
each iteration. In the next section some of the parameters of equation (22) will be eval- 
uated and a final form useful for computation will be presented. 


Obtaining Working Equations for Analog 

Before equation (1) can be solved in the analog computer, it is necessary to consider 
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scaling the problem and calibrating the analog. It is convenient to discuss calibration 
briefly at this point in connection with the evaluation of various parameters occurring in 
equation (22). 

Consider the x-component of equation (1) written in terms of the real system 



dt 
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or converting to the model-size thrustor, 
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Define 
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where 6 is a small potential difference measured between two closely spaced points 

y III 

in the electrolyte along the x-axis and K c is a constant. Actually, $ x m is the signal 
input to the analog computer measured in the electrolyte by the two x-component sam- 
pling probes. Equation (1) then becomes 



= r)aK c 6 


x, m 


( 26 ) 


If the time scaling ratio 0 = = dr^/dt^ is used, equation (26) may be written as 
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( 27 ) 


A: 7«K c 6 xm 

^ e 2 

The relation between the distance variable x and the computer output in volts from 
equation (6) is 


x= — = — (28) 

T F FS 

where X is the distance variable expressed in volts and S is the plotting scale factor 
in volts per meter. Therefore, 




(29) 


All the components of the right side of equation (29) are constant for any given problem 
except 5 x m , so that 


d?X = B6 


dr? 


x, m 


(30) 


where 


B = 


7)FSaK c 


(31) 


The significance of B in equation (30) is that it represents the overall gain of the signal 
6 x m that is fed to the com P uter - The value of B determines the plotting speed and is 
adjusted by a potentiometer until this speed is optimized (i. e. , plotting too fast will 
create waves in the electrolyte and cause distortion; plotting too slow will allow small 
errors to multiply in the electronic integrators). 

The y-components of the equations of motion are developed in exactly the same man- 
ner as the x-component. Therefore, 
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In order to evaluate I in equation (22), it is 
necessary to obtain 0. From equation (31), 


/r?FS«K c 
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(33) 


and from equation (25), 
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dVr 


3x f 6 
i x, m 


(34) 


In the calibration dV^/dx^ = V c /x c where V c is the voltage across the calibration elec- 
trodes (see fig. 12), x„ is the distance between them, and 6 * 6 is the signal fed 

to the computer from the sampling probes on the surface of the electrolyte between the 
electrodes, so that K c may be evaluated from 

K c =— (35) 

X n 6 n 

c c 

The values of the terms in equation (35) are the same for both x- and y-components. 
Finally, by substituting for /3 in equation (22) and by using equation (35), the current to 
be injected in a simulation volume of the model system is given by 


1 = 


4V2cB 1 / 2 Txj/ 2 e c 1 / 2 V?/ 2 

wyv/2 



(36) 


Equation (36) shows that the space-charge simulation current is independent of the scale 
of the models since F is contained only in B, which is an operationally determined 
constant. Once a converged solution is obtained for the trajectories, the current density 
for the solution may be calculated by using the plane diode approximation, by choosing 
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Solution 



Figure 13. - Comparison of space-charge- limited-flow poten- 
tial distributions between coaxial cylinders under different 
polarization effects. 


the proper value for rj, and by substituting the thrustor voltages and dimensions de- 
sired. The generality of equation (36) is obvious since it applies to all two-dimensional 
configurations. Equation (36) will be referred to in the next section where some example 
solutions are discussed. 


ILLUSTRATIVE EXAMPLES 
Initial Test of Analog and Discussion of Solutions 

A cylindrical diode model for which the analytic space-charge-limited flow results 
are known (ref. 10) was used to check the quality of the results obtained by the 
electrolytic -tank analog. Aluminum electrodes were used initially in this test because 
of their availability. They were then gold-plated to check the polarization effects. Com- 
parison of the analytic with the analog space-charge-limited potential distributions is 
shown in figure 13 for the two different surface materials. Because some aluminum was 
exposed, the electrodes were very difficult to keep clean for more than a few minutes in 
the water. 

For the gold-plated electrodes, which were the best, the maximum error in trajec- 
tory position for the test case was 3 percent where the error is defined as the trajectory 
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deviation from the analytic solution at its endpoint divided by the total trajectory length. 

Some of the solutions to be presented in this section will be compared with solutions 
obtained from a digital computer program (ref. 9). The computer furnishes more de- 
tailed results for a final analysis. As will be shown in this section, though, the electro- 
lytic tank analog is more versatile in initial studies. Once trends are established and 
the optimum electrode configuration is approached by using the electrolytic -tank analog, 
the more accurate digital program may be desired for the final stages of solutions. 

The two models to be discussed in this section are possible configurations for a 
contact- ionization type of electrostatic thrustor. These configurations are based on the 
idealized design concepts given in reference 11. 


Divergent-Flow Model 

As discussed in reference 11, the divergent-flow concept has the promise of high 
current densities at the emitter and reduced densities in the area of critical charge ex- 
change (aperture area). A model of the divergent-flow electrostatic thrustor in the elec- 
trolytic tank is shown in figure 4 (p. 5). The space-charge-free solution is shown in 
figure 14(a). This solution is a starting point toward a practical design for the thrustor. 
Although figure 14(a) shows only a space-charge-free solution of the model, it serves as 
a reference for studying configuration variations (e. g. , the trend of the magnitude of 
the current density can be observed by noting the changes in position of an equipotential 
line near the emitting surfaces of various models). These variations are of interest for 
the purpose of improving the average current density possible from the thrustor and 
focusing the beam more satisfactorily. A flat electrode downstream of the accelerator 
aperture is maintained at effectively zero potential (i. e. , ground reference) to simulate 
an approximate exhaust beam neutralization plane. This in turn allows simulation of 
accel-decel operation of the model. In the configuration shown, the ratio used was 
^A^net = ^ where q> ^ is the difference between the emitter and the accelerator and 
</? ne t * s net potential difference between the emitter and the ground. 

An important factor that must be minimized in the design of a thrustor is charge- 
exchange erosion of the accelerator (ref. 1). Charge exchange may occur when a fast 
moving ion collides with a slow neutral particle. After the collision the slow particle 
has acquired the charge and, depending on its originating position, may strike the elec- 
trodes, causing erosion and reduced thrustor lifetime. As discussed in reference 2, the 
electrolytic tank analog is well suited for studying the charge -exchange ion trajectories 
and patterns of electrode erosion. Although charge-exchange studies made in a Laplace 
potential field do not directly apply to the final solution of a model, they may be used in 
the manner previously discussed to predict trends. Figures 14(b) to (d) show space- 
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(b) Configuration 2. 

Figure 14. - Divergent -flow configurations of space-charge-free solution in electrolytic tank. Focuser voltage in tank, 0 volt; accelerator 
voltage in tank, 6 volts; neutralization plane voltage in tank, 2 volts. 
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(a) Configuration 1 of space -charge-free solution. Average current density 
at ion-emitting surface, 159 amperes per square meter. 




(b) Configuration 1 of space -charge -limited solution. Average current density at 
ion -emitting surface, 128 amperes per square meter. 

Figure 16. - Circular -flow configurations of various solutions in electrolytic tank. 
Focuser simulated voltage, 2.8 kilovolts; accelerator simulated voltage, -7.2 kilo- 
volts; neutralization plane simulated voltage, 0 kilovolt. 


charge-free solutions for some con- 
figuration variations of the original 
model. The dashed lines in fig- 
ures 14(b) to (d) represent charge- 
exchange ion trajectories. The 
small circles are typical points 
where the charge exchanges are as- 
sumed to occur. The charge - 
exchange ions formed at these points 
start out with zero initial velocity. 

To obtain these trajectories, the 
field- sampling probe is positioned at 
the desired point and the trajectories 
are begun with no initial condition 
voltage (i. e. , zero initial velocity). 
Many of these trajectories can be 
started at various points in the beam 
envelope so that a boundary may be 
established for the region where det- 
rimental charge exchanges can oc- 
cur. If the neutral efflux, primary 
ion density, and electric potential at 
the starting point are known, these 
trajectories provide a means for de- 
termining erosion patterns on the 
thrustor electrode and calculating 
electrode life expectancy. 

The space-charge-limited-flow 
solution for the model shown in fig- 
ure 14(a) was obtained by using the 
electrolytic tank and also the digital 
computer program of reference 9. 
For the electrolytic -tank solution 
(fig. 15(a)), the emitter was divided 
into eight equal parts to establish the 
start of the stream tubes. The tra- 
jectories are started in the center of 
each stream tube so that the trajec- 
tory nearest the accelerator is not 
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idi Configuration 3 of space charge-free solution. Average current density at ion- 
emitting surface, 175 amperes per square meter. 

Figure 16. - Concluded. Circular-flow configurations of various solutions in elec- 
trolytic tank. Focuser simulated voltage, 2.8 kilovolts; accelerator simulated 
voltage, -7.2 kilovolts; neutralization plane simulated voltage, 0 kilovolt. 


the beam boundary. The short 
dashed line represents this boundary. 
Figure 15(b) shows the digital solu- 
tion for the model. In each case an 
accel-decel ratio (4>^/<l> net ) °* ^ was 
used. A comparison between the two 
solutions can be made based on the 
amount of the aperture occupied by 
the beam envelope and the average 
emitter current density predicted. 
The two solutions differ in beam 
envelope size by 7. 3 percent at the 
accelerator aperture. The lines 
traced in figure 15(b) represent 
stream-tube boundaries, while those 
in figure 15(a) represent the center- 
lines of the stream tubes. The pre- 
dicted average current density at the 
emitter j differs by 9. 3 percent. 
The difference in predicted current 
density is probably due to inaccura- 
cies in the measurement of £ f - 
(eq. (21)) from the electrolytic -tank 
trajectory plot. 


Circular-Flow Model 

As discussed in reference 11, 
the theoretical advantage of the 
circular -flow configuration is that 
the electrodes provide thermal 
shielding, thus reducing heat losses 
(fig. 16). With this configuration 
the hot ionizer does not have a direct 
line to space. 

As in the previous case, the 
Laplace solution shown in fig- 
ure 16(a) is the starting point for 
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analysis of the concept. Figure 16(b) shows the space-charge-limited-flow solution for 
this configuration. The space charge reduced the average current density of the emitter 
by 19. 5 percent. This reduction is quite significant at high current densities. Also, 
figure 16(b) shows that the ion beam has been spread and shifted toward the elliptical ac- 
celerator electrode. The beam boundary barely misses the accelerator electrode. This 
close approach is undesirable because mechanical variations will probably result in im- 
pingement in the real thrustor. The ratios of the voltages on the various electrodes for 
these cases are those discussed in reference 11. 

The effect on the space-charge-free solution of displacing the circular arc electrode 
nearest the beam is shown in figure 16(c). As can be seen, this provides more clearance 
between the beam envelope and the electrodes. Variations that simplify the design of the 
thrustor are also of interest. As an example, figure 16(d) shows a variation in which 
the large, curved outer electrode was replaced by a smaller, simpler one. The inner 
curved electrode was replaced to its original position. Also, the neutralization plane or 
simulated virtual ground was removed to observe the effect. The beam envelope pro- 
duced by this variation is narrower than the previous configuration. The average pre- 
dicted current density at the emitter increased 10 percent over the original model. 

On several occasions during the solutions for the two thrustor geometries discussed 
in this section, an iteration on the solution for both geometries was made in the same 
day. This illustrates the convenience of the removable-pin tray design and the versatil- 
ity of this tank analog as an electrostatic -thrustor design tool. 


CONCLUDING REMARKS 

The Lewis Research Center electrolytic tank analog has been used to analyze 
electrostatic -thrustor geometries that present a particularly wide variety of two- 
dimensional space-charge-flow problems. The uniform spacing of the current- injection 
pins in this analog was found to provide a sufficiently accurate space-charge simulation 
for these configurations. This uniform spacing and the use of easily removable pins 
greatly reduced the time required to change models and make the high standard of clean- 
liness required for accurate solutions easier to maintain. The plane diode approxima- 
tion method used for calculating the current flow in the stream tubes permitted generali- 
zation of the space-charge injection-current calculations, thereby further increasing the 
versatility of the analog. The circuitry developed for differencing, multiplying, and 
filtering the computer input signals increased the reliability and simplified the electronic 
requirements of the analog. 

In the initial investigation of a problem, the space -charge-free solutions were found 
to be extremely useful guides to point out design deficiencies and establish trends for 
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reaching final electrode configurations and spacings in a space-charge field. More de- 
tailed solutions for the final stages of analysis are possible with digital computer pro- 
grams. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 26, 1965. 


APPENDIX A 


SYMBOLS 


fraction of stream-tube sec- 
tion within simulation vol- 
ume 

amplitudes of sinusoidal 
function 

emitting surface area, sq m 
constants 


dv 

element of volume 

X,Y 

coordinates, V 

x, y, z 

Cartesian coordinates 

a 

ratio of field intensities 

0 

time scaling factor 

6 

small potential difference, V 


square wave distortion fac- 
tor, V 

electric field intensity, V/m 

dimensional scale factor 
space-charge simulation 
current, A 

current, A 

current density, A/sq m 

calibration constant 

linear distance 

mass of ion, kg 

number of time increments 

square wave amplitude, V 

starting point of trajectory 

electric charge, C 

charge of ion, C 

plotter speed, V/m 

element of surface 

time per increment, 
sec/space 

real time, sec 
electric potential, V 


e Q permittivity of free space 

77 charge- mass ratio of ion, C/kg 

8 phase difference 

p space-charge density, C/cu m 

a conductivity of electrolyte, 

mho/m 

r scaled time, sec 

total accelerating potential, V 

4>net ne t accelerating potential, V 

u> frequency, cps 

Subscripts: 
c calibration 

e emitter 

f model-size thrustor 

i i in stream tube 

m analog model 

r real model 

v simulation volume 

0 initial value at starting point of 

trajectory (appendix B) 
Superscript: 

C ) average 



APPENDIX B 


DETAILS OF INJECTION-PIN-CURRENT CALCULATIONS 

Equation (36) is used to calculate the pin currents for space-charge simulation. The 
value of for each trajectory stream tube is measured by plotting an equipotential line 
near the emitter (e. g. , see fig. 10, p. 13). The conductivity of the electrolyte is meas- 
ured immediately before the pin-current calculations are made. 

The first approximation to the pin-current value is too high because of the method 
employed to find the space-charge-limited flow (ref. 9). If the currents are decreased 
too much, a converged solution may require many additional iterations; if not enough, 
either a negative gradient is set up (an unrealistic condition) or the succeeding iterations 
will oscillate (i. e. , trajectories will converge, diverge, etc.). 

The following calculations are made from the sample plot shown in figure 17. The 
first number of simulation- volume time increments (from emitter to trajectory starting 
point P ) in square A2 are determined by the quotient of the time it took the particle to 
travel from the emitter to point P r divided by the time T between timing marks in 
seconds. Some typical values for the terms in equation (36) are the following: 


y 
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Conductivity of electrolyte, a, mho/m 2. 0X10“ ^ 

Constant B 0.2 

Time per increment, T, sec/space 0. 5 

Distance between calibration electrodes, x , m 0.1016 

Small potential difference between calibration electrodes, 6 , V 4. 0 

Electric potential of model-size thrustor, V^, V 0. 1 

Emitter surface area of model-size thrustor, a^, sq m 3.68x10“^ 

Calibration electric potential, V , V 4. 0 

Plotter speed, S, m/V 197. 0 


so that equation (36) becomes 


I = 7. 3X10' 10 J]N V> i\q 2 i (Bl) 

i 

Tables I to III illustrate the calculation of the injection current needed for the two 
pins shown in figure 17. Thus from table II and equation (Bl), the pin currents for this 
sample iteration are those given in table HI. The voltage dropped across the 10-kilohm 
series resistor must therefore be 3. 51 volts for square A2 and 0. 758 volt for square B2. 


TABLE I. - STARTING POINTS, INITIAL VELOCITIES, AND C 2 . 

1 


Stream 

tube, 

i 

Coordi- 
nate of 
starting 
point, 

x o 

Initial 
velocity 
in x- 

direction, 

x 0 

Time 

lapse 

to 

starting 

point, 

t 0’ 

sec 

Coordi- 
nate of 
starting 
point, 

^o 

Initial 
velocity 
in y- 

direction, 

yo 

Time 

lapse 

to 

starting 

point, 

t 0 » 

sec 

Plane 
diode 
of sepa- 
ration, 

h, i> 

m 

r 2 
f, i’ 

-2 

m 

5 

0. 0077 

-0. 0072 

2.2 

0. 0011 

-0. 0009 

2.2 

10. 0X10" 3 

G 

6 

.0076 

-.0070 

2.2 

. 0014 

-.0011 

2.3 

9.0 

1 

7 

. 0076 

-.0066 

2.3 

. 0016 

-.0012 

2.4 

8.0 

KM 

8 

.0075 

-.0060 

2. 5 

. 0019 

-.0013 

2.6 

7.0 

2. 04 

9 

. 0073 

-.0053 

2. 8 

. 0022 

-. 0013 

2.8 

6. 5 

2.37 
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TABLE H. - CALCULATION OfVn .A. ft: 2 . 

V, 1 11,1 

i 


Square 

Stream 

tube, 

i 

Number of 
time incre- 
ments in i** 1 
stream tube 
for volume 

v, 

N v,i 

Prorating 

factor, 

A i 

N v, i A i 

Wm 

i 

A2 

5 

(2. 2/T) +5.9 

0. 15 

1.21 

1.21X10 4 

48. 15X10 4 


6 

(2. 2/T) + 6. 1 

.98 

8. 13 

10.00 




(2. 3/T) + 6.3 

1.00 

8.60 

13.42 




(2. 5/T) + 6. 6 

. 90 

8. 19 

16.70 




2.8 + 6.8 

.30 

2.88 

6.82 


B2 

5 

2.6 

0.60 

1.56 

1. 56x1 0 4 

10. 38X10 4 


6 

2.6 

1. 00 

2.60 

3.20 




2.8 

. 96 

2.69 

4.20 




2.8 

.25 

.70 

1.42 



TABLE HI. - INJECTION 
PIN CURRENTS FOR 
TWO SAMPLES 


Square 

Pin current, 


amp 

A2 

3. 51X10" 4 

B2 

7. 58X10" 5 



















APPENDIX C 


EVALUATION OF VARIOUS ELECTRODE MATERIALS 


T 


1 


Figure 18. - Square wave obtained by use of 
brass electrodes. 


The model electrode material must be nonreactive to the tap water. For this rea- 
son, brass and stainless steel were considered for the material. Brass is easily ma- 
chined and has the desired electrical properties; however, it was found that with brass 
electrodes, distortion of the electric field occurred because of polarization at the inter- 
face of the electrode and the electrolyte. The extent of distortion of this kind can be 
measured by applying a square wave signal across the electrodes and reading the signal 

on an oscilloscope. The leading edge of the 
square wave shown in figure 18 indicates 
distortion. The relative distortion of var- 
ious electrode materials and surface condi- 
tions, such as contamination caused by 
cleaning with alcohol or the presence of oil 
or fingerprints, can be compared by obtain- 
ing for each material and surface condition 
the ratio d/2P. 

The test was conducted by using plane 
diode electrode configurations with 8 inches 
of spacing, and a 5. 4- volt square wave was 
applied. The signal is measured with a 
probe immersed in the electrolyte between 
the electrodes at various distances from the 
grounded electrode. The distortion d/2P 
is the percent reduction in the area under 
the square wave. The data from the com- 
parison is shown in figure 19. It can be seen 
that best results were obtained with clean 
stainless steel. It is also noteworthy, how- 
ever, that the gold-plated brass gave good results that were also very insensitive to 
surface cleanliness. These electrodes have advantages, but the gold-plating is time con- 
suming and costly. Stainless steel is more difficult to fabricate but requires no surface 
preparation other than keeping the critical surfaces of the electrodes free from oils and 
fingerprints. 



Normalized distance from grounded 
electrode 

Figure 19. - Polarization effects using various electrode 
materials. 
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